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Abstract. We extend our previous approach (Eur. Phys. J. B, 74, 63(2010)) to modeling correlated elec- 
tronic states and the metal-insulator transition by applying the so-called statistically consistent Gutzwiller 
approximation (SGA) to carry out self-consistent calculations of the renormalized single-particle Wannier 
ryj ■ functions in the correlated state. The transition to the Mott-Hubbard insulating state at temperature 

T = is of weak first order even if antiferromagnetism is disregarded. The magnitude of the introduced 
d self-consistent magnetic correlation field is calculated and shown to lead to a small magnetic moment in 

the magnetically uniform state. Realistic value of the applied magnetic field has a minor influence on the 
metallic-state characteristics near the Mott-Hubbard lcalization threshold. The whole analysis has been 
' O , carried out for an extended Hubbard model on a simple cubic (SC) lattice and the evolution of physical 

Ch ■ properties is analyzed as a function of the lattice parameter for the renormalized ls-type Wannier functions. 

Quantum critical scaling of the selected physical properties is analyzed as a function of the lattice constant 
i^i , R — s> R c = 4.1ao, where R c is the critical value for metal-insulator transition and ao = 0.53A is the Bohr 

radius. A critical pressure for metallization of solid atomic hydrogen is estimated and is ~ 10 2 GPa. 

PACS. 71.30.+h Insulator-metal transitions - 71.27.+a Strongly correlated electron systems - 71.10.Fd 
Lattice fermion models - 62.50.-p High-pressure effects in solids and liquids 

CO 

' 1 Motivation electron states and the theory of fermionic quantum liq- 

CN| , uids with delocalized states, for which (quasi)momentum 

O ■ One of the important quantum-mechanical problems in P = hk characterizes quasiparticle states. A principal dy- 

■ both the solid-state [fll] and the optical-lattice g] sys- namic quantity driving the transition is the particle den- 

tems is the localization-delocalization transitions of fer- sit y (or interatomic distance for hxed number of particles, 

mionic states which is called the Mott or Mott-Hubbard as i s the case here). 

transition. In electron systems it corresponds to the delo- A second impetus to the physics of these phenomena 

calization of atomic states (usually of 3d or 4/ type) and has been provided by the introduction by Anderson [5] and 

S3 ' the formation of a Fermi liquid of heavy quasiparticles Hubbard [7] of second-quantization or quantum-particle 

composed mainly of the transforming electrons, irrespec- language, with the help of which the revised Mott-Hubbard 

tively of the fact that there may be other valence elec- transition can be analyzed in terms of microscopic param- 

trons present also there [3]. In the extreme version, the eters. In the simplest, half-filled single-band model the rel- 

transition is driven solely by the interparticle interaction, evant microscopic parameter is the ratio U/W, where U 

in which the presence of lattice plays only a secondary is the magnitude of the intraatomic (Hubbard) interac- 

role. In that situation, the localization-delocalization tran- tion and W = 2z\t\ is the bare bandwidth (i.e., that for 

sition is called the Wigner transition [5]. The description the uncorrected particles), with z being the coordination 

of these two transitions minimally provides a bridge be- number (i.e., the number of nearest neighbors) and \t\ is 

tween the atomic physics with the localized single-particle the magnitude of intersite transfer (hopping) of individ- 

ual fermions. It is amazing that a similar type of approach 

a e-mail: kadzielawa@th.if.uj.edu.pl can be formulated for both fermions and the bosons, in the 

b corresponding author: ufspalek@if.uj.edu.pl latter situation in the optical-lattice situation [3J. 
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The principal question is how to combine the Mott [TJ 
and the Hubbard [8] aspects of this quantum phase tran- 
sition in a purely electronic system. Once achived, the 
whole description can be analyzed as a function of the 
lattice parameter (or interatomic distance. R). In the se- 
ries of papers [MTU] we have formulated such an approach 
starting from the Gutzwiller-ansatz approximation (GA) 
for the extended Hubbard model, with a simultaneous 
readjustment of the single-particle Wannier functions in 
the correlated state. The method provides, among oth- 
ers, a direct connection of the Mott criterion for local- 

1/3 

ization/delocalization (n c as ~ 0.2) with that of Hub- 
bard (U « W). As an extra bonus coming from such 
a formulation we obtain the quantum-critical behavior 
of single-particle wave- function size a -1 [TO], as well as 
the detailed evolution of the correlated metallic state into 
the Mott-Hubbard insulator. One principal methodolog- 
ical advantage of the present formulation is that, in dis- 
tinction to the LDA+U PUDl] or LDA+DMFT [T3j ap- 
proaches our formulation avoids the problem of double 
counting of electron-electron interaction. However, unlike 
LDA+U or LDA+DMFT methods, the present approach 
is still on the stage of modeling only the simplest (one- 
band) systems. Nonetheless, it may provide a formally 
proper starting point for more complex situations such 
as many-band systems. Our formulations can certainly 
be also reformulated for Bose-Hubbard optical-lattice sys- 
tems. It is unique in the sense of discussing of quantum 
critical behavior of single-particle wave-function charac- 
teristics such as the inverse wave-function size (the inverse 
effective Bohr radius). 

As a concrete application of our results we calculate 
the critical pressure for the metallization of solid atomic 
hydrogen with the electronic correlations included within 
our renormalized mean-field theory. 

The structure of the paper is as follows. In Sec. [5] we 
characterize briefly the method used in the paper and the 
modification of our previous approach (S][TUJ. In Sec. [3] 
we analyze in detail the numerical results obtained with 
the help of the so-called Statistically Consistent Gutzwiller 
Approximation (SGA). As a physical application, we also 
provide there the estimate of the critical pressure for the 
solid atomic hydrogen metallization. Section0]contains an 
outlook with summary of main results and a brief discus- 
sion of possible extensions. 



2 Model and Method Applied 
2.1 Starting Hamiltonian 

We start with the Extended Hubbard Hamiltonian for Is 
hydrogenic-like system 

H=6 a J2 n *+ Y t i3 a ur a j< r + U Y n it n ii 

' **" _ 1 (1) 

i<j i<j 



where tij is the hopping integral, U the intraatomic in- 
teraction magnitude, e a the atomic energy per site, and 
Vi on -i on corresponds to classical Coulomb interaction be- 
tween two H + ions, equal (in atomic units) to 



Vi 



ton— ion \' j 



(fj - n) 



(2) 



By following [JJ], we introduce N e = ^ rii - the total 
number of electrons, and Srii = m — 1 as the deviation 
from neutral-atom configuration. We can now rearrange 
the interatomic interaction in the following manner 



i<3 i<3 



Y K ij n i n j =y^ j K i j{n i - l)(nj - 1) 

i<j 

i<j i<j 

i<j 



(3) 



N 

i<j 

By introducing now effective atomic energy per site in 
the form 



^ 4 

we can rewrite Hamiltonian JTJ in the following way 

U = e a ff Y H i + Y ^J a L a ja + U Y n A n ^ 

+ KijSnMj. 



(4) 



(5) 



We add — a o~\giiBH a ni a - a simple magnetic Zee- 
mann term , where g is the Lande factor, \i& the Bohr 
magneton and H a the external magnetic field. By intro- 
ducing the reduced magnetic field h = ^g^BH ai we obtain 
our starting Hamiltonian 



H =e a // Y n% + Y <'J"^"," + U Y n ^ n ^ 

i i^j,cr i 



-Y K i3 Sn i$ n 3 ~ h Y' 



(0) 



This lattice Hamiltonian describing the system of ls- 
type states in a solid contains the following microscopic 
parameters: t^ , U, Kij and the band filling n. Addition- 
ally, close to the metal-insulator boundary we can assume 
that (SniSnj) ~ 0. 



2.2 Incorporation of single-particle wave function 
optimization 

The microscopic parameters of ((SJ) can be expressed via 
the single-particle Wannier functions in the following man- 
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Uj = {wi\Hi \wj) 



U = ( WiW. 



Kij = ( WiWj 



\ri - r 2 \ 

Q 2 



ri - r 2 \ 
(w i \H 1 \uii) . 



WiW, 



WiWj 



(7) 



where Hi is the Hamiltonian for a single particle in 
the system, and e 2 /|ri — r 2 \ is interparticle interaction. 
The numerical value is obtained by approximating first 
the Wannier functions Wi = Wi (r) = Wi (r — Ri) by Is 
Slater orbitals and those, in turn, by series of the Gaussian 
functions, i.e., 



w 



(r)=^(r)- 7 ^^(r), 

3=1 



&i(r) 



a e - a \ r -R t \ 
7T 



(8) 



3 

a 2 



EM? 



-r?\r—Ri\ 



j3 and 7 = \/\ — j3 2 parameters depend explicitly on 
the overlap integrals of !^ functions and z is the num- 
ber of nearest neighbors. Parameters B a and r a are de- 
rived by minimizing energy of single atom (Hamiltonian 
H a = — v 2 — 2|r — Ri\ ). p is the number of Gaussian 
functions used. Parameter a is found by minimizing the 
system energy of the trial correlated state (see belows). 
The difference with our previous approach QJ5] is that we 
include the statistical consistency conditions, as discussed 
next. 



2.3 Statistically-Consistent Gutzwiller Approximation 
(SGA) 

To obtain optimal value of the inverse size a given in- 
traatomic distance R we have to obtain the system en- 
ergy. It was shown |15| that the Gutzwiller Approximation 
(GA) does not always provide the variational results con- 
sistent with those obtained from the corresponding self- 
consistent equations. To assure this consistency we mini- 
mize the GA free-energy functional J- supplemented with 
two additional molecular fields A m and A n , coupled with 
m and n respectively: 



T 



(SGA) _ 



(9) 



+ A (X n n + X m m + Ud 2 + fin) 
where the trial eigenvalues E^ ^ are: 



A„ 



= q a £k - <j(h + A m ) - X n ) , 



where d 2 = {n^n^) , and 

q a = (y/ (n a - d 2 ) (1 - n a - n w + d 2 ) 
+ dyjn-a- d 2 ^j jn„ (1 - n a ) 



(11) 



is the band narrowing renormalization factor and £k 
is dispersion relation for bare particles (here taken for 
simple cubic structure, = It (cos k x + cos k y + cos k z )). 
The eigenvalues Ek„ are obtained by Fourier transform of 
the effective GA Hamiltonian 



n = e a ff n ™ + t i3^ a L a ja+ 

ia ija 



(12) 



additionally supplemented with the Lagrange-multiplier 
constrains 

- X m (ntj - m) - A n J~] (rij - n) , (13) 

i i 

where m ?; = - n^, m = (m,), m = n it + n^, and 
n = (rii}. 

The operator /C = H — ^2n(X m mi + X n ni) + A(X m m + 
X n n) plays the role of the effective Hamiltonian, in which 
the mean fields (m, d 2 ) and the Lagrange multipliers (A m , 
A n ), as well as fx, are all determined variationally, in addi- 
tion to the wave function parameter a. In the next section 
we analyze in detail the physical results obtained by min- 
imizing though that procedure the functional 



2.4 Overview of numerical methods 

Numerical analysis was carried out with the help of method 
different from that used in [MTU] by introducing the lower- 
level minimization for each single-particle basis optimiza- 
tion step. In other words - for each and every step of min- 
imizing energy with respect to the reverse function size a, 
where Golden Section Search |16| was empirically proven 
to be the most efficient, there is introduced a new mini- 
mization of functional J^( SGA ) with respect to the double 
occupancy d, the magnetization m, the chemical potential 
jj, and the molecular fields X m and A„. The latter pro- 
cedure was carried using GSL - GNU Scientific Library 



Qhttp : //www . gnu . org/ sof t ware/gsl/) , with the order 
of magnitude of zero-precision 10 _s and all the following 
up calculations with the double precision. 

Due to new minimization step and new parameter (ex- 
ternal magnetic field H a ), the numerical complexity in- 
creases by two orders of magnitude, and it enforces new 
optimization. To decrease the computing time we chose 
the basis of three Gaussians (p = 3 in ©) instead of seven 
(as in [9ll0| ). After comparing the results for constant ex- 
ternal magnetic field H a for p = 3 and p = 7 we observe 
no qualitative change of behavior. Below we discuss the 
basic physical properties obtained within SGA and just 




Fig. 1. Ground state energy (per atom) of the metallic state 
(PM, x) for R < R c and the insulating (PI, +) for R > R c , as 
a function of interatomic distance R (in units of Bohr radius 
ao). Inset: detailed representation of the first-order PM — ► PI 
transition near R — R c ~ 4.1ao. The upper curve for R < R c 
represents the energy of the unstable PI state. Note that as 
Eg > —lRy, the lattice can only be stabilized by the external 
pressure (see sec. 13. 3p . 



discussed numerical procedure, as well as compare them 
with those obtained previously without the statistical con- 
sistency pnu]. 



Fig. 2. Double occupancy probability d 2 = (nifn^) versus R. 
Note a weak discontinuous jump to zero at R = R c ~ 4.1ao, 
as compared to the continuous evolution for PM — > PI of the 
Gutzwiller approximation obtained previously [9lll0|. 



> 



3 Results and discussion 




3.1 Ground-state properties 



Fig. 3. Microscopic parameters: hopping integral |i| and the 
Hubbard interaction parameter U, both as a function of R. In- 
set: U/W ratio vs R. The U/W ratio for R = R c is (U/W) c ~ 
1.18 . The vertical dashed line marks the PM — > PI tran- 
sition point. In the large- R limit U reaches the atomic value 
Uat = (5/4) Ry. 



The calculated ground state energy Eg as a function of in- 
teratomic distance R (lattice parameter) is shown in Fig. 
[TJ In the inset a detailed dependence of Eq in the tran- 
sition regime R « R c — 4.1ao is displayed. For the sake 
of comparison, the energy of PI state (for H a = 0) has 
been also shown for R < R c , where this phase is not sta- 
ble. The principal difference with out previous Gutzwiller 
ansatz (GA) results is that in the present (SGA) approach 
the transition is weakly discontinuous, as one can see ex- 
plicitly from the circumstance that (IEq /dR, which is pro- 
portional to the internal pressure, is discontinuous. 

In Fig. [2] we plot the double occupancy probability 
d 2 = (riifriii) versus R and again see a weak discontinuity 
at R = R c . The present SGA results are compared with 
the previous GA results of Spalek et al. [TU] . 




3.4 3.6 3.8 4 4.2 4.4 



R(a ) 

Fig. 4. Correlation-induced effective magnetic field A m vs. R 
for selected values of applied magnetic field H a - For H a — 0, 
A m vanishes at PM — PI boundary (see inset for details). The 
presence of the applied field triggers A m 7^ in the insulating 
state. 
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The independent evolution of the microscopic param- 
eters are shown in Figs. [3] and [4] In Fig. [3] we plot the val- 
ues of U and the nearest- neighbor hopping magnitude \t\, 
whereas in the inset the U/W ratio, with W = 2z\t\ being 
the bare bandwidth, is shown in the regime R « R c . Note 
that at R c = 4.1eto the U /W is of the order of unity, in ac- 
cordance with the results obtained when U jW is treated 
as a free parameter [5]. The atomic limit is reached effec- 
tively when \t\ sw 0, i.e. for R ~ 6.5ao. The value of U 
in that limit is U = (5/4) Ry. Amazingly, one should also 

note that at the transition R^a^ 1 = ric Ob ~ 0.25, in 
accordance with the original criterion due to Mott [2] for 
the Mott localization (ric = 1/Rc is the particle density 
at R = R c and as = or 1 in the effective Bohr radius at 
that point) . This approach does allow to relate directly the 
Mott and the Hubbard criteria for localization of fermions. 
This was achieved by readjusting the Wannier functions, 
determining t and U parameters in the correlated state. 

In Fig. 5] we exhibit the evolution of the consistency 
field A m , which plays the role of the self-consistent (cor- 
relation) field, also for selected nonzero applied magnetic 
field. One should note that the lowermost curve corre- 
sponds to the case with H a = 0. The nonzero value of the 
effective field in the metallic phase will introduce a small 
but nonzero value of the spin magnetic moment in the 
uniform (weakly magnetic) phase as discussed below. The 
situation is highlighted in the inset to this figure. How- 
ever, one should emphasize that the field A m is dependent 
on the phase discussed and in the antiferromagnetic state 
it takes the form of a staggered field [T7]. The latter phase 
will not be analyzed in detail here, as the first and the 
foremost aim of this paper is to underline the first-order 
nature of the M — > I transition if the statistical consis- 
tency conditions are properly taken into account, as well 
as the quantum scaling of the single-particle atomic wave 
function inverse size ov _1 . The corresponding behavior of 
the renormalized but less pronounced (cf. [TO]), so we will 
not discuss it in detail here. 

To visualize directly the type of magnetism accompa- 
nying the metal- insulator transition (MIT) , we have plot- 
ted in Fig. [5] the zero-field magnetic susceptibility versus 
R (the inset illustrates an overall behavior) . The principal 
feature of \ quasi-divergence for R R c is roughly inde- 
pendent of the applied field magnitude specified. It does 
not follow exactly the Brinkman-Rice dependence [TS] as 
here the unrenormalized (xo) value is strongly R depen- 
dent as well, so x cannot be parametrized by U/W ratio 
only. For the sake of completeness we plot in Fig. |5] the 
magnetization curve at the localization threshold. The de- 
pendence roughly emulates the Brillouin curve for the lo- 
calized particles of spin 1/2. However, one has to take 
into account a nonlinear increase of the molecular field for 
nonzero H a value, (cf. Fig. [4j) . In the inset to this figure 
we provide the uniform magnetic-moment magnitude as 
a function of R if a uniform A m ^ is assumed. In the 
metallic state the spontaneous value of spin polarization 
is very small and practically field independent. 

For the sake of completeness, we also provide in Tab.[T] 
the numerical values of the selected quantities vs. R. Note 



100 



if 50 




Fig. 5. Differential static magnetic susceptibility vs. R for se- 
lected values of the applied magnetic field. The \ divergence at 
R c — 4.1ao accompanies the PM — PI transition and is associ- 
ated with localization of the itinerant electrons when R — » R c . 
Overall \ behavior in the metallic state does not depend much 
on the value of H a - Inset: double logarithmic plot x(-R) showing 
absence of any simple exponential type of scaling. 




50 100 150 

applied field H a (T) 



200 



Fig. 6. Magnetic moment m — (n^ — n^) as a function of 
the applied magnetic field at the critical interatomic distance 
R c = 4.1ao. Inset: a residual uniform moment in the metallic 
state vs. R near R c , induced by the correlation field when it is 
assumed as spatially uniform. 



that we compared the values of Eq in SGA and GA. The 
slightly higher values of Eg obtained in the present (SGA) 
case should not be surprising for R away from R c . It is re- 
warding that they are lower as R c is approaching, since 
then the results are more realistic, i.e. the tight-binding 
approximation works better. The last column describes 
the effective mass enhancement (or 1 = due to the 
interparticle correlations. One should note that the bare 
band value tub is also dependent on R. Additionally, for 
the half-filled band case (and only then) the mass enhance- 
ment is not explicitly spin-direction (a — ±1) dependent. 
The scaling with R of the selected just discussed quantities 
is carried out next. 
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Table 1. Computed quantities in SGA (except additional results for ground-state energy added, calculated in GA), all as a 
function of the lattice parameter, for the simple cubic structure. For their definition, see main text. The values, if not specified 
explicitly, are in rydbergs (Ry). The x(0) value for R ^ R c is infinite. 



R(an) 


Eg GA 




t 


u 


cy 1 ( G,r\\ 


d 1 


A 






3.25 


-0.8640 


-0.8644 


-0.2409 


1.4996 


0.9474 


0.152774 


0.015884 


0.1809 


1.17728 


3.50 


-0.8814 


-0.8816 


-0.1773 


1.4749 


0.9220 


0.120128 


0.012641 


2.0598 


1.36818 


4.00 


-0.9136 


-0.9136 


-0.1098 


1.4152 


0.9200 


0.048886 


0.006084 


22.6577 


2.82781 


4.05 


-0.9171 


-0.9171 


-0.1046 


1.4139 


0.9175 


0.038973 


0.004256 


47.3562 


3.47235 


4.09 


-0.9200 




-0.1005 


1.4140 


0.9147 


0.030193 


0.027281 


253.7567 


4.40375 


4.10 


-0.9209 


-0.9207 


-0.0995 


1.4143 


0.9138 


0.000000 


0.000000 


oo 


oo 


4.20 


-0.9315 


-0.9288 


-0.0896 


1.4217 


0.9021 


0.000000 


0.000000 






4.50 


-0.9544 


-0.9517 


-0.0705 


1.3742 


0.9263 


0.000000 


0.000000 






5.00 


-0.9760 


-0.9732 


-0.0471 


1.3200 


0.9556 


0.000000 


0.000000 






7.00 


-0.9939 




-0.0082 


1.2504 


0.9972 


0.000000 


0.000000 






oo 


-1.0000 


-1.0000 


0.0000 


1.2500 


1.0000 


0.000000 


0.000000 







3.2 Scaling laws of physical parameters 

The question arises whether the scaling properties pro- 
posed earlier [TU] of selected physical quantities near the 
quantum critical point (located at R = R c ), where \ 1S 
divergent, take place also within the present SGA ap- 
proach. To address this question, we have plotted in Fig. 
[7] the relative ground state energy per site \{Eg{R) — 
E G (R C ))/E G (R C )\ = \(E - E c )/E c \ vs. relative inter- 
atomic distance (R c — R)/R c for R — > R c — 0. As before, 
we observe an almost linear scaling \(E — Ec)/Ec\ 
= A((R C - i?)/i? c ) 7 , since 7 = 1.035 ± 0.009 and ad- 
ditionally, A = 0.66 ± 0.04. One should note that this 
type of scaling for the ground state energy near QCP 
is very different from the corresponding behavior in the 
atomic limit R 3> i? c , where the corresponding energy 
per site is roughly dominated by the Coulomb contribu- 
tion, so Eg ~ R . This type of scaling means that the 
first derivative is nonzero and constant for R — > R c . This 
feature will be taken upon and interpreted in Section l3~51 



10 

1 

0.1 
0.01 



LU 
LU 



0.001 - 

0.0001 r 

1e-05 



H .0349+0.0089, A=0.658+0.036 

data obtained ■ 




1e-05 0.0001 0.001 0.01 
(R c -R)/R c 



0.1 



Fig. 7. Scaling of the relative ground state energy near R — 
R c . For detailed explanation see main text. 
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1 |r 

0.1 
0.01 
0.001 
0.0001 



1e-05 



y=1 .021 +0.044, A=0.45±0.13 
data obtained 



1e-05 0.0001 



0.001 0.01 
(R c -R)/R, 



0.1 



Fig. 8. Scaling of the inverse orbital size a near R = R c . For 
detailed explanation see main text. 



In Fig. [5] we plot the relative inverse renormalized 
atomic-orbital size | (a(R)—a(R c ))/a(Rc) \ = \(a — a c )/a c \ 
vs. (R c — R)/R c . Again, we have an almost linear scaling, 
though the behavior is not as regular as before [TU]. It 
is thus more legible to display the actual behavior of the 
renormalized orbital size a(R) vs. R, as shown in Fig. [5J 
The quantum critical behavior corresponds to the cusp- 
like feature around i? c , where we observe a clear discon- 
tinuity of da 1 1 dR at R = R c . The most important feature 
of this figure is the fact that we observe this type of quan- 
tum critical behavior for the renormalized inverse size of 
the single-particle wave function. We are not aware of any 
other result of that type appearing in the literature, apart 
from the analogical result obtained by us earlier within the 
GA [TDj. Note also, as already said earlier, this type of crit- 
ical behavior will translate into the same type of behavior 
for the renormalized Wannier functions. It is important to 
emphasize that this quantum critical behavior of the prin- 
cipal characteristic (a) of the renormalized wave function 
is independent of the fact whether the actual two-particle 
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characteristic - d 2 vanishes continuously (as in GA) or 
discontinuously (as in SGA), as well as whether the states 
are paramagnetic (PM, PI) or have a small "parasitic" 
magnetic moment. This analysis should certainly be ex- 
tended and carried out for the antiferromagnetic (AFM, 
AFI) cases, but that requires a separate treatment, as it 
will certain a large number of additional variational pa- 
rameters. This is certainly planned as the next stage of 
work. 



1.12 



1.08 



1.06 




1.04 - 



1.02 - 



Fig. 9. Inverse atomic-orbital size vs. R. The critical behavior 
near R = R c = 4.1ao. The dependence is practically indepen- 
dent of the applied field. Note that the cusp has its maximum 
slightly above R c . 



3.47 



2.83 



2.82 
1.84 



1.83 




R=3.8 a, 



o 



50 



Ha (T) 



100 



150 



Fig. 10. Effective mass enhancement as a function of mag- 
netic field for the situation with one electron per atom. In this 
case it is spin direction independent and dependent only on 
the interatomic distance. Note that this enhancement does not 
include the change of the bare mass which increases also with 
the increasing R. 



Finally, we would like to discuss at this point two im- 
portant physical characteristics of our correlated single 



narrow-band system. First of them is exhibited in Fig. 
[TOland shows the effective mass enhancement q~ x vs. H a 
for selected values of R close to R c , but on the metal- 
lic side. We observe the two important features. First, 
the quasiparticle mass is spin-independent in this half- 
filled band case, unlike in the non-half- filled case [TTjll20) . 
Also, the mass enhancement due to the correlations in- 
creases in SGA relatively slowly compared to that in the 
parametrized GA approach [2T] as R — > R c — 0. Second, 
this slow increase of q~ 1 when confronted with the rather 
fast dependence of the susceptibility (cf. Fig. [SJ Tab. [TJ) 
means that the magnetic contribution to the renormalized 
susceptibility (corresponding to the "renormalized Stoner 
part", cf. [18,21]) has an essential contribution to the sus- 
ceptibility. Also, as mentioned earlier, there may be an es- 
sential contribution due to the bare bandwidth W narrow- 
ing with the increasing R. In effect, the resultant behavior 
of quantities such as x or the linear specific-heat coefficient 
7 = 7o9 _1 in the strong-correlation limit (here U ps W) is 
much more subtle than it was discussed within the pure 
parametrized model (GA) picture, where the renormaliza- 
tion of the Wannier functions, as well as the correlation 
fields A m and A„, are absent. Simply put, our approach 
goes beyond the parametrized-model approach as it pro- 
vides the correlated system evolution as a function of ex- 
perimentally accessible parameter - the lattice parameter. 




Fig. 11. Hubbard gap evolution E gap — U — W as a function 
of interatomic distance. Note a continuous dependence near R c 
even though the transition is weakly discontinuous (cf. Figs.[T] 
[2]). As the gap is increasing rapidly with the increasing R (R> 
Rc), the Mott insulators fir R > 4.5ao cannot be regarded as 
wide-gap semiconductors. 



Finally, in Fig.[TT]we plot the value of the Mott-Hubbard 



gap E g 



U — W versus R. Obviously, in the metallic 



phase it does not appear (for the sake of completeness, 
we show it as having a negative value) . Above the transi- 
tion it reaches relatively fast with the increasing R > R c 
the asymptotic value E gap ~ U. A weak cusp-like behav- 
ior should be noted just above R = R c . This feature is 
reminiscent of the cusp-like behavior of the renormalized 
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Waimier functions. It would be interesting to study the 
pressure dependence of the gap for the Mott insulator lo- 
cated just above the transition point and determine its 
scaling with (R — R c ) / R c . It does not scale linearly, as the 
quantities earlier. 



3.3 Supplement: Application to solid atomic hydrogen: 
estimated of the critical pressure for metallization 

Our formulation represents a fully microscopic model of 
the solid atomic hydrogen undergoing a Mott insulator- 
metal transition. This model differs from the standard 
treatment |22j . in which the phase diagram is treated as 
a function of microscopic parameter U/W, since we in- 
clude here a fully self-consistent procedure of evaluating 
the renormalized-by-correlations wave functions. In effect, 
we can calculate explicitly the critical pressure for metal- 
lization. In this approach the external pressure is the fac- 
tor stabilizing the system in a particular phase (M or /). 
This task can be carried out by using the classical defi- 
nition of pressure p as the force F, applied to make the 
present atomic solid stable, over the area A. This force 
is obtained by differentiating F = — \7rEg\- The corre- 
sponding external pressure we have to exert on the system 
in order to stabilize the crystal as a function of interatomic 
distance (with energy Eq per site and the elementary cell 
area A/N = R 2 ) is plotted in Fig. 1131 Note that the phys- 
ical meaning has the critical pressure for metallization as 
the point R = R c is the terminal point of stability of 
the almost localized Fermi liquid. It differs from other re- 
sults |23] as those results concern a molecular (H2) solid; 
for detailed comment see e.g. |24| . Note that the pres- 
sure behavior p(R) for R — > R c traces the difference on 
the slopes dE/dR for metallic (R < R c ) and insulating 
(R > R c ) regimes. Roughly, it shows a critical (divergent) 
behavior at R c . To determine that more accurately we re- 
quire a refined numerical analysis. Nonetheless, the trend 
is clear and provides a promising starting point for a fur- 
ther analysis. 

One must emphasize that this discussion is based on 
an implicit assumption that the positive ions (protons in 
this case) form a lattice in both the metallic (M) and 
insulating (/) states. This may turn only a rather crude 
approximation on the possible role of their zero-point mo- 
tion in increasing the value of the critical pressure. Addi- 
tionally, a possible transition to solid-liquid plasma as an 
alternative scenario, should certainly not be ruled out (for 
the relevant discussion of that point e.g. [2"51l26] ). 



4 Outlook 

We have provided ground-state principal characteristics of 
the correlated narrow-band system involving first-principle 
variational calculations of the renormalized Wannier func- 
tions of ls-type in the correlated state described by the 
Statistically Consistent Gutzwiller Approximation (SGA). 
Our model is fully microscopic as it contains an explicit 
calculations of microscopic parameters of the extended 



Q_ 

o 



ext. pressure p=|-V R E G |/R 
Insulator 




7.5 8 



Fig. 12. External pressure (in GPa = 10 9 Pa) one has to exert 
in order to stabilize the crystal vs. interatomic distance (with 
cell area A/N = R 2 ). Note two critical values of pressure: 
pc, 1 = 97.7GPa required in the insulating state and pc,M = 
62.6GPa in the metallic state. The lattice becomes very rigid 
as R —¥ R c - The dotted line marks naive extrapolation M —¥ I. 



Hubbard model. In this manner, we can track the cor- 
related Fermi-liquid evolution with the increasing lattice 
parameter R into the Mott-Hubbard insulator. As an ex- 
tra bonus, we obtain an estimate of the critical pressure 
for the metallization of solid atomic hydrogen which is 
roughly equal to lOOGPa. Obviously, the solid atomic 
hydrogen is a molecular crystal up to much higher pres- 
sures [2"2"ll2"5] , but it may be possible that the states of this 
type can be synthesized in the optical-lattice experiments. 

We would like to emphasize that our model calcula- 
tions contain two attractive features from a general method- 
ological point of view of the computational solid state 
physics. First, they do not include twice the Coulomb in- 
teractions among electrons, as it is the case for LDA+U or 
LDA+DMFT methods. Those latter methods have their 
own merits as they can be (and are) applied to com- 
plex real materials. But the correctness of the present ap- 
proach is still worth mentioning, even when applied so far 
to the mode system only. Second, within the present ap- 
proach the calculation of the renormalized characteristics 
of the single-particle wave function and of the interelec- 
tronic correlations, are both treated on the same footing. 
This is the crucial feature of the Mott-Hubbard systems, 
for which the single-particle energy (e.g. the renormalized 
band energy) is of the same magnitude as the interparticle 
(Coulomb in this case) interaction. The approach of the 
present type should be extended to the multiband situa- 
tion to discuss the realistic strongly correlated materials 
(the magnetic oxides such as V2O3) evolving, when vary- 
ing the lattice parameter (i.e., applying the pressure), and 
not only as a function of the microscopic parameters such 
as U/W, as they vary very rapidly in the vicinity of the 
Mott-Hubbard transition. However, this program execu- 
tion may represent a long road ahead. 

One has to mention that the present approach requires 
some other basic extensions. First of all, as any Gutzwiller- 
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type approach, it does not include explicitly the intersite 
kinetic exchange appearing deeply in the Mott insulat- 
ing state (for R substantially larger than R c ). The trace 
of this interaction is coded in the correlation fields A m , 
as has been shown in a related context before [T7]- Sec- 
ond, the case with partial band filling (n < 1) should be 
also treated. Third and most importantly, the antiferro- 
magnetic state should be included in the analysis. This 
is our plan for the near future. However, we believe that 
this simple analysis shows up in a clear manner the quan- 
tum critical behavior of the renormalized-by-correlations 
single-particle wave function, not obscured by the compli- 
cated magnetic structure (i.e., that with a staggered mag- 
netic moment). Fourth and finally, the approach, even in 
the present single-band case, should be extended to the 
close packed lattices such as fee and hep, as this are the 
typical structures for metals. But then, one has to include 
at least first two hopping integrals (between the nearest 
and the next- nearest neighbors). Such treatment is cer- 
tainly tractable. Nonetheless, we believe that our first es- 
timate of the critical pressure for the metallization of solid 
atomic hydrogen, carried out for the sc structure, provides 
a promising starting point. 

Finally, one can extend our analysis to nonzero tem- 
peratures in a straightforward manner, starting from the 
free-energy functional ^ . This should be carried out sep- 
arately and the results can be compared with GA results 
providing the first-order line ending in a classical critical 
point, as well as the reentrant metallic behavior [19 20 27 . 
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